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We report measurements of conditional Eulerian and Lagrangian structure functions in order to 
assess the effects of non-universal properties of the large scales on the small scales in turbulence. We 
study a Im x Im x 1.5m flow between oscillating grids which produces Rx = 285 while containing 
regions of nearly homogeneous and highly inhomogeneous turbulence. Large data sets of three- 
dimensional tracer particle velocities have been collected using stereoscopic high speed cameras 
with real-time image compression technology. Eulerian and Lagrangian structure functions are 
measured in both homogeneous and inhomogeneous regions of the flow. We condition the structure 
functions on the instantaneous large scale velocity or on the grid phase. At all scales, the structure 
functions depend strongly on the large scale velocity, but are independent of the grid phase. We 
see clear signatures of inhomogeneity near the oscillating grids, but even in the homogeneous region 
in the center we see a surprisingly strong dependence on the large scale velocity that remains 
at all scales. Previous work has shown that similar correlations extend to very high Reynolds 
numbers. Comprehensive measurements of these effects in a laboratory flow provide a powerful tool 
for assessing the effects of shear, inhomogeneity and intermittency of the large scales on the small 
scales in turbulence. 



I. INTRODUCTION 

Many of the most powerful insights in the study of 
fluid turbulence are rooted in the idea of an energy cas- 
cade where the chaotic process of transferring energy to 
smaller scales allows the small scales to become universal 
and independent of the details of the forcing mechanism 
at the large scales. However, careful examination of many 
small scale statistics in different flows has shown that 
the reality of turbulence is quite a bit more complicated. 
Some statistics such as the scaling exponents of Eulerian 
and Lagrangian struc ture functions are nearly identical 
in different flow^i'^ISl g^^; other small scale statistics 
such as the coefficients in scaling lawPl^ or the scalar 
derivative skewnesJ^ show dependence on the properties 
of the large scales up to the largest Reynolds numbers 
measured. 

A traditional approach to deal with dependence on the 
large scales has been to classify flows, and allow that 
there might be differences between categories of flows 
such as free shear flows (jets, mixing layers, etc), wall 
bounded shear flows (boundary layers, channel flows, 
etc), or isotropic turbulence (wind tunnel grid turbulence 
or numerical simulations in a box with periodic bound- 
ary conditions, etc). A careful empirical comparison of 
statistics between different flows could then show which 
properties of the small scales are truly independent of the 
large scales. 

Recent developments in experimental tools have al- 
lowed for the measurement of small scale Lagrangian 
statistic^. In contrast to flows designed for Eulerian 
measurements such as hot-wire anemometry, Lagrangian 
measurements do not require large mean velocity in or- 
der to use Taylor's frozen flow hypothesis. In fact the 
opposite is desired since small mean velocity allows a 



particle to be tracked in an observation volume for the 
longest possible time. This has led to significantly dif- 
ferent flow designs for Lagrangian measurements. Many 
new flows have been introduced which have complex large 
scales and are difficult to place into traditional categories 
which were originally meant for flows designed for hot- 
wire anemometry. Two widely used examples of flows 
with low mean velo city are counter rotating disk d^ l ^ * ^ -^ 
and oscillating grid^^nHEll^ yg^ newer generation of 
flows are currently under study including the random 
jet arra}E^, corner stirred tankl^^, Lagrangian exploration 
module (LEM), and radial acoustic jets which each have 
unique large scales. 

Initial work on Lagrangian measurements in flows with 
complex large scales and a small mean velocity has pri- 
marily assumed that the small scale statistics of inter- 
est are independent of the large scale forcing of the flow. 
This assumption has been tested in several cases by thor- 
ough quantitative comparison of small scale statistics in 
different flows. Comparison of acceleration probability 
distribution functions (pdf) in direct numerical simula- 
tions (DNS) and counter rota ting di sk experiments are 
found to be flow independent I 9 | i6 | i7| ^ jj-^ addition, the 
scaling exponents of the Lagrangian structure functions 
have been compared between DNS and experiment, and 
found to be in close agreement!^. Much more work is 
needed to determine the degree to which the large scales 
of different flows affect various Lagrangian statistics. 

However, there is a more direct way to evaluate 
whether the small scales in a flow are independent of the 
large scales: the small scale measurements can be condi- 
tioned on a measurement of the state of the large scales. 
Two previous studies have shown strong dependence of 
the small scales of Eulerian structure functions on the 
instantaneous velocity, which is dominated by the large 
scales. Praskovsky et alP^ extensively study interactions 
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between the large scales and inertial range scales in two 
high Reynolds number wind tunnel flows. Strong correla- 
tions between the large scales and the velocity structure 
functions are found at all length scales. They interpret 
this as being consistent with the correct application of 
Kolmogorov theory with a fluctuating energy injection 
at the large scales. Sreenivasan and Dhruval^ measured 
Eulerian velocity structure functions from atmospheric 
boundary layer data for R\ > 10"', some of the largest 
Reynolds numbers ever measured. They find that the 
structure functions conditioned on the large scale veloc- 
ity show a strong dependence, and they show that DNS 
and grid turbulence measurements show almost no de- 
pendence. They attribute the dependence to large scale 
shear, and show how to remove the effect to improve 
power law scaling. There is evidence that other small 
scale statistics show conditional dependence on the large 
scales. The acceleration variance shows a strong depen- 
de nce w hen it is conditioned on the large scale veloc- 

One challenge in discussing interactions between large 
scales and small scales is the very non-universal nature 
of the large scales. Each flow has a unique set of large 
scales, which may depend on time, geometry, or driving 
parameters. So it has been difficult to isolate the aspects 
of the large scale ffow that are affecting the small scales. 
Anisotropy is the aspect that is best understood. Ex- 
tensive work has identified persistent anisotropy at small 
scales even at very high Reynolds number&fi^, and anal- 
ysis using spherical tenso r deco mposition has placed this 
problem on solid footing mMl. However, this is not the 
only effect of the large scales. Here we wish to distinguish 
two additional aspects of the large scales that are partic- 
ularly important. Inhomogeneity is the spatial variation 
of statistics. Large scale intermittency is temporal ffuc- 
tuations on time scales longer than the eddy turnover 
time, L/u. Both inhomogeneity and large scale intermit- 
tency often occur together in real ffows, but are distinct 
properties since ffows can be conceived that have each 
without the other. For example, a homogeneous turbu- 
lent ffow in DNS can have large scale intermittency by 
having the energy injection varied in time. 

In this paper we present a comprehensive set of mea- 
surements of the dependence of Eulerian and Lagrangian 
velocity structure functions conditioned on the large scale 
velocity. We use a ffow between two oscillating grids 
which is relatively homogeneous in the central region, but 
has large inhomogeneity near the grids. This allows us to 
isolate the signatures of different properties of the large 
scales. We ffnd clear signatures of inhomogeneity, but a 
significant part of the dependence of the structure func- 
tions on the large scale velocity seems to be the result of 
large scale intermittency. A better understanding of this 
dependence on non-universal large scales will help in the 
identification of universal statistics, and the comparison 
of different ffows. 



II. EXPERIMENT 

This work is based on optically tracking passive tracer 
particles seeded in a turbulent ffow agitated by two oscil- 
lating grids as shown in Fig. [T] For clear measurements of 
1-D inhomogeneity a large system is needed to create sig- 
nificant separation for locally homogeneous and inhomo- 
geneous regions, and to create sufficiently high Reynolds 
numbers. In order to study large scale effects conditional 
statistics were analyzed, which required large data sets 
(> 10^ particle pairs). Storage, speed, and budget con- 
cerns led to the development of real time image compres- 
sion circuits^. These devices enabled nearly endless data 
acquisition for a nominal cost. 

A. Experimental Apparatus 

Turbulence was generated between two identical octag- 
onal grids oscillating in phase in an octagonal Plexiglas 
tank that is Im x Im x 1.5m and filled with approxi- 
mately 1,100 liters (300 gallons) of filtered, degassed wa- 
ter. The grids have 8 cm mesh size, 36% solidity, and 
were evenly spaced from the top and bottom of the tank 
with a 56.2 cm spacing between grids and a 1 cm gap 
between grid and tank walls. The stroke was 12 cm 
peak to peak, powered by an 11 kW motor. A typical 
grid frequency was 3 Hz, but was raised to 5 Hz to in- 
vestigate Reynolds number dependence. Water cooling 
maintains the temperature at ± O.IC during each run. 
Neutrally buoyant 136 /xm diameter polystyrene tracer 
particles were added to the flow until approximately 50 
were seen by each camera. This particle density was cho- 
sen to maximize data per frame while minimizing track- 
ing errors. Particle density could greatly increase with 
the planned addition of two more cameras^® One dif- 
ficulty of oscillating grid fiows is that vibrations from the 
oscillatory drive can couple to the camera supports and 
degrade imaging accuracy. We mounted a custom cam- 
era support on an optical table to minimize vibrations. 
Air bubble suppression was an additional concern. We 
developed a method to keep all water seals and bearings 
sufficiently wet to maintain an air tight seal. 

B. Detection 

These data were acquired using 3D PTV (particle 
tracking velocimetry) measurements using two Bassler 
A504K video cameras capable of 1280 x 1024 pixel res- 
olution at 500 frames per second (a data rate of approxi- 
mately 625 MB per second per camera). Recording such 
a high data rate is a significant technological hurdle. A 
typical system would store data in 4 GB of video RAM, 
so that one run would last just 7 seconds before wait- 
ing approximately 7 minutes for the data to download 
to hard disk. We use an image compression circuit to 
threshold images in real-time so that only pixels above 
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FIG. 1: Experimental apparatus diagram. Two oscillating 
grids were held 56.2 cm apart in a 1,100 liter octagonal prism 
Plexiglas tank. Two high speed cameras are used to stereo- 
scopically image chosen regions of the tank in order to record 
3D particle positions. Illumination is provided by a Nd:YAG 
laser with 50 W average power. 



a user defined brightness limit are regarded as particle 
data and retained while the dark background pixels are 
discarded^. This technique produces a dynamic data 
compression factor of 100-1000, which enables continu- 
ous data collection and storage to hard disk. 

Our first implementation of the image compression cir- 
cuits have faced two major challenges. First, the simple 
thresholding compression reduces particle center accu- 
racy. However, particle finding accuracy is typically de- 
graded by only 0.1 pixel, which is typically less than the 
uncertainty in particle finding from unthresholded im- 
ages. Secondly, because frame number information was 
created and recorded separately on each computer, any 
operating system delay can lead to frames lost and tim- 
ing mismatch between the cameras. For the measure- 
ments in this paper, frame number errors were corrected 
in post processing. Updated versions of the image com- 
pression circuit have solved this problem by including 
camera frame number in the data stream the computers 
record. 

Particles are illuminated using a 532 nm pulsed 
Nd:YAG laser with 50 W average power. The beam 
was expanded to create an illumination volume approx- 
imately 7cm X 4cm x 5cm. Images were processed to 
find the center of each particle as seen by each camera 
and then stereomatched to find the 3D position in real 
space. Stereomatching was accurate to approximately 
11 fim (0.08 particle diameters or 0.2 pixels). At this level 
of accuracy it is essential to have a very good calibration 
of camera position parameters to use for stereomatching. 
We start with a traditional calibration to obtain initial 
3D stereomatching'^. We then use known stereomatched 



pairs from the two cameras, and run a non- linear opti- 
mization to minimize the stereomatching error and find 
optimal camera position parameter^2ll 



III. RESULTS 
A. Characterizing the Flow 

We define a characteristic velocity by u = {{uiUi) /3)-^/^ 
and a characteristic length scale by L = v?/e where 
e is the energy dissipation rate per unit mass defined 
in section IIIC For the center region u — cm/s, 

8.3 cm/s. 



L = 9.0 cm, and for the near grid region u 
L — 4.5 cm. The Taylor Reynolds number, R\ = 
{IbuL/vY^^, (where v is the kinematic viscosity) ranges 
from 285 for 3 Hz grid frequency to 380 for 5 Hz grid 
frequency in the center. Near the grid at 3Hz R\ = 230. 
The Kolmogorov length and time scales are rj — 140/im, 
T,, = 20 ms in the center region, and 77 = 94 /im, — 8.8 
ms in the near grid region. 

Figure [2]^ a) shows the mean vertical velocity as a func- 
tion of the vertical position along the central axis of the 
tank. The top and bottom grids are separated by 56.2 
cm, approximately 7 L. In Fig. [2] the dot-dashed line in- 
dicates the maximum amplitude of the bottom grid, 22.1 
cm below the center of the tank. Data was collected at 
5 separate heights in order to measure the complete flow 
profile from the center of the tank to the bottom grid. 
Mapping the bottom half of the tank is sufficient because 
the geometrical symmetry produces a mirror image above 
the midplane. The two volumes which we will focus on 
throughout this paper are bounded by the dashed lines, 
and will be referred to as the center (C), and near grid 
(NG) observation volumes. At this grid separation dis- 
tance the mean flow traces four torii, two above and two 
below the center plane of the tank, as shown in the sketch 
in figure [3] (drawn to scale). In the large central region, 
the effect of the mean flow is to pump highly energetic 
fluid from the region near the grid towards the center of 
the tank. In figure |2ja), there are two points where the 
mean vertical velocity reaches zero: one near the center, 
the other 18 cm from the center just below the near grid 
observation volume. The existence of this second stag- 
nation point and reverse circulation region depicted in 
Fig. [3] is a common feature in mean flows generated by 
oscillations^^. In the following measurements, the small 
mean velocity has been subtracted so that we study the 
fluctuating velocity field. 

Figure [2j^b) shows the vertical velocity variance along 
the central axis as a function of vertical position. The 
velocity variance is large near the grid and quickly falls 
off towards the center where it is nearly homogeneous. 
The center and near grid observation volumes were cho- 
sen to provide a contrast between the large homogeneous 
region in the center and the much more inhomogeneous 
region near the grid. In the center, the variance of the 
velocity is homogeneous for several L in either direction. 
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FIG. 2: Mean and variance of vertical velocity along the 
central vertical axis of the tank. Grid frequency 3 Hz, grid 
separation distance 56.2 cm. The dot-dash line represents the 
grid height at majcimum amplitude. The remainder of the 
paper will focus on measurements in two regions designated 
by the vertical dashed lines: one at the center of the tank and 
one near the grid. 



The velocity variance ranges moderately in the near grid 
observation volume, and enormously within one L below 
this region. In Fig. [2] deviations from a smooth curve 
are not due to statistical uncertainty, but are a result of 
patching 5 calibrated regions together with the majority 
of error coming from measuring absolute position in the 
tank. 

It is interesting to note that we made measurements in 
a flow with a smaller grid separation of 35 cm and found 
that the Reynolds number in the center was lower. The 
characteristic velocity in the center did increase due to 
the closer proximity of the grids, but L was reduced by a 
larger amount resulting in approximately 8% decrease in 
R\. The reason for the unexpected decrease in Reynolds 
number is a reversal of the mean velocity compared with 
larger grid separations. For larger grid separation dis- 
tances, energetic fluid from near the grids is carried to 
the center by the mean flow. However, at 35 cm grid 
separation the mean velocity reverses which results in a 
lower Reynolds number in the center. 



B. Structure Functions 

To measure Eulerian structure functions we first find 
the instantaneous longitudinal velocity difference be- 
tween two particles a distance r apart Aur = [u(x) — 




FIG. 3: Scale diagram of 56cm x 100cm area between grids 
showing the mean circulation torii which are nearly rotation- 
ally symmetric about the central vertical axis. Center and 
Near Grid observation volumes are drawn in dashed lines, 
which shows the relative size and position of the observation 
volumes in figure |2] Horizontal dot-dashed lines represent 
range of motion of the top and bottom grids. 



u(x -I- r)]i, where the L subscript denotes the longitu- 
dinal component, found by projecting the 3D velocity 
difference vector onto the vector connecting the two par- 
ticles. The longitudinal structure functions are defined 
as Dp — {{Aur)P) where p represents the order of the 
structure function and the brackets represent the ensem- 
ble average. In the inertial range, Kolmogorov (1941) 
gives 



(A<)=Cf)(er)f/3, 



(1) 



where the Cp are Eulerian Kolmogorov constants and 
e is the energy dissipation rate. 

Figures |4] and [5] show the measured second and third 
order longitudinal velocity structure functions with the 
straight thin lines representing Kolmogorov's prediction 
from Eq. [T| The insets show the structure functions com- 
pensated by Eq. [T] At Rx = 285, any scahng range is 
very limited, but the plateaus can be used to estimate 
the inertial range. 

Lagrangian structure functions were measured from 
temporal velocity differences along a particle trajectory. 
The velocity difference now becomes Au^. = u(t) — u{t + 
r), where r is the time interval between measurements. 
We use the vertical component of the velocity for La- 
grangian velocity differences throughout this paper, al- 
though results for the other components are similar. For 
Lagrangian structure functions Kolmogorov (1941) pre- 
dicts 



(A<) = Cf)(er)P/2 



(2) 
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FIG. 4: Eulerian second order longitudinal velocity structure 
function shown as a function of pair separation r scaled by the 
Kolmogorov length rj. The inset shows this data compensated 
by Eq.[l]for p=2. 



FIG. 5: Eulerian third order longitudinal velocity structure 
function. The inset shows this data compensated by Eq. [T]for 
p=3. 



C. Energy Dissipation Rate Measurement 

The energy dissipation rate e is an important value 
throughout this analysis, it is worth a moment to discuss 
how it is determined. Limitations in particle density pre- 
clude direct measurement via the definition 

e = 2v{sijSij) 

with 

1 / dui duj \ 
'■^ 2 \dxj dxi ) ' 

Instead we utilize Kolmogorov 's 4/5 law: Eq. M with 
p = 3 where the coefficient C3 — —4/5 can be de- 
rived from the Navier-Stokes equations. We identify the 
inertial range with the plateau in the compensated third 
order structure function (Fig. p] inset). The inertial range 
is chosen to be 25 to 91 r/Ty (735 to 1.3 cm). If the same 
inertial range is used in the second order structure func- 
tion, the energy dissipation rate determined from it (us- 
ing the empirical coefficient Cj^-* = 2.0)'^ is within 3% 
of the value calculated from the third order. 



D. Phase Dependence 

A simple energy cascade has constant energy input at 
the largest length scales. An obvious departure from con- 
stant energy input is the oscillating grid driving mecha- 
nism. The sinusoidal motion of the grid directly corre- 
sponds to energy input with periodic time dependence. 
It seems likely that such a strongly periodic energy in- 
put would have a signature throughout the whole energy 
cascade. 

The method we employ throughout this work to detect 
signatures of the large scales is to condition various statis- 
tics on some measurement of the state of the large scales. 
In this case, we condition on the phase of the grid mo- 
tion, (p. Conditioning instantaneous single particle statis- 
tics such as the mean and variance of the velocity shows 
some sinusoidal dependence on grid phase. For exam- 
ple, in the center the conditional variance, ((u— (u))^|(/)), 
varies by 1% over the cycle of the grid. In the near grid 
region, the conditional variance varies by 10%. The mean 
vertical velocity in the center, {u\(j)), varies by 0.8 cm/s 
over a cycle of the grid which is 10% of the standard 
deviation. Near the grid, the conditional mean velocity 
varies by 2 cm/s, which is 20% of the standard deviation 
at that location. 

Figure [6] shows the compensated second order Ion- 



6 



gitudinal structure functions conditioned on the phase, 
(Au^|0). In the center of the flow (Fig. |6ja)) the struc- 
ture functions have essentially no change with phase. 
Near the grid (Fig. ^h)) there is a slight dependence 
on grid phase. To emphasize the differences between 
structure functions at different phases, we compensated 
the structure functions by a single energy dissipation 
rate in each figure, £3 = 24.6 cm^/s^ in the center and 
£3 = 131 cm^/s'^ near the grid. These values were de- 
termined when the grid is in mid amplitude (the third 
bin) . The good collapse of the structure functions at all 
phases across the entire range of r shows the minimal 
dependence of the small scales on the large scale peri- 
odicity of the flow created by the oscillating grids. One 
possible source of dependence of small scale statistics on 
large scales has been shown to be minimal. 
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FIG. 6: (Color Online) Second order compensated veloc- 
ity structure functions conditioned on grid phase. The col- 
lapse shows the very weak phase dependence, a) center of the 
tank (C). b) near the grid (NG). Zero and 2n phase represents 
grid at lowest possible amplitude, (p: + = - 2tt/5, * = 2tv/5 
- 47r/5, o = 47r/5 - 6tv/5, A = 6tv/5 - Stv/5 , □ = 8tt/5 - 2tv. 



E. Dependence on Large Scale Velocity 

1. Eulerian structure functions conditioned on the large 
scale velocity: center region 

A more revealing dependence on the large scales of the 
flow is found by conditioning the velocity structure func- 
tions directly on the large scale velocity. A convenient 
measurable quantity that reflects the local instantaneous 
state of the large scales is the average velocity of the 
particle pair used for the structure function, defined as 
T,Uz = {uz{x) + Uz{x + r))/2. Alternatively, conditioning 
on the average velocity of many particles, not just one 
pair, was studied and found to have very similar results, 
but we choose to focus on Su^ because it can be more 
easily measured and does not depend on the observa- 
tion volume and seeding density. Additional conditioning 
quantities will be discussed in section |III G[ 

Figure [7]^a) shows the second order Eulerian velocity 
structure function conditioned on Eu^. The smallest val- 
ues of the structure function correspond to pair velocities 
near zero, represented by o, while large results in 

larger values of the structure functions. For the bins we 
chose, the structure function conditioned on large values 
of TiUz is nearly twice the value when conditioned on Sit^ 
near zero. 

Figure ^h) shows the data in Figure [Tf^a) compen- 
sated by Kolmogorov inertial range scaling. The func- 
tional forms are quite similar, confirming the impression 
from Fig. that all length scales are affected simi- 
larly by the instantaneous state of the large scales. In 
Fig. [rj^b) we used a different energy dissipation rate, e„^ 
to compensate each of the 5 individual large scale verti- 
cal velocity bins. This insures all conditions plateau at 
approximately the same value, and allows for direct com- 
parison of the functional forms of the conditional struc- 
ture functions. 

The strong dependence of the conditional structure 
functions on the large scale velocity at all scales reveals 
that the small scales are not statistically independent of 
the large scales in this flow. There is not any detectable 
trend toward the smaller scales becoming less dependent 
on the large scale velocity than somewhat larger scales. 



2. Lagrangian structure functions conditioned on the large 
scale velocity: center region 

In much the same way we can evaluate the conditional 
Lagrangian structure functions. Figure |8ja) shows the 
second order Lagrangian structure function conditioned 
on the vertical component of the large scale velocity, Eu^ . 
Here we find TtUz by averaging the velocity of the particle 
at the two times used to determine Aur- The conditional 
structure functions for different large scale velocity are 
different by a factor of about 2.5, and they remain nearly 
parallel throughout the entire time range. 
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FIG. 8: (Color online) Second order Lagrangian velocity 
structure function conditioned on instantaneous velocity (ver- 
tical component) in the center of the tank, a) Uncompensated 
structure functions, b) Individually compensated to have the 
peak values match. Symbols represent the following dimen- 
sionless vertical velocities, Eii^/ -^Z {u|): + = 3.1 to 1.9, * = 
1.9 to 0.62, o = 0.62 to -0.62, A = -0.62 to -1.9 , □ = -1.9 to 
-3.1. 



FIG. 7: (Color online) Second order velocity structure func- 
tion conditioned on particle pair velocity (vertical component) 
in the center of the tank, a) Uncompensated structure func- 
tions b) Individually compensated by the energy dissipation 
rate for each conditional data set. Symbols represent the fol- 
lowing dimensionless vertical velocities, Sit^/ {u|): -f — 4.2 
to 2.5, * = 2.5 to 0.84, o = 0.84 to -0.84, A = -0.84 to -2.5, 
□ = -2.5 to -4.2. 



Figure |8[b) shows the second order conditional La- 
grangian structure function compensated by Eq. [2] where 
e is individually chosen so the maxima of all of the con- 
ditioned structure functions coincide. This aids compar- 
ison of the functional forms of the conditioned structure 
functions. Again, the functional form is nearly identi- 
cal for different large scale velocities, indicating that the 
large scales affect all time scales in the same way. There 
may be a small trend towards larger values of the com- 
pensated Lagrangian structure functions at small times 
when the magnitude of the large scale velocity is large. 

It should be noted that there is a bias present in La- 
grangian measurements that is not present in Eulerian 
measurements. A sample of measured Lagrangian tra- 



jectories is biased towards low velocity particles since 
the high velocity particles are more likely to have left 
the measurement volume. This bias becomes larger for 
larger r. Berg et al..^'^ have studied this bias and find 
that it can be quite large for typical experimental con- 
ditions. We quantified this bias in our data by measur- 
ing the Lagrangian structure functions using trajectories 
that remained inside artificially restricted measurement 
volumes. From a simple extrapolation of the dependence 
on the size of the artificial detection volume, we esti- 
mate that our experimental Lagrangian structure func- 
tions underestimate the true value by 17% for r = 8r^. 
This is roughly consistent with the size of the error we 
expect based on the critical time lag defined in Ref.'SSl 
Note that we have not performed the compensation they 
recommend and we are roughly translating their uncom- 
pensated results. Because of this bias, we will focus at- 



tention on T < 10t„. As we'll discuss in section IIIF3 



the dependence of the conditioned Lagrangian structure 
functions on the large scale velocity does not seem to be 
significantly influenced by this bias. 
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3. Eulerian structure functions conditioned on the large 
scale velocity: near grid region 



By comparing separate regions of the tank we are able 
to explore the effects of inhomogeneity on this condi- 
tional dependence. Figure [9] shows the Eulerian struc- 
ture functions, similar to Fig. [7] but with data collected 
in the inhomogeneous region near the bottom grid (NG) . 
The separation between Eulerian structure function con- 
ditions has doubled to approximately a factor of four. 
Note the different ordering of the structure functions. 
The up-down symmetry is now broken. Fluid travel- 
ing upwards (* symbols) has a large structure function 
while fluid traveling downward with the same magnitude 
of vertical velocity (A symbols) has the lowest value of 
the structure function. We interpret this as highly en- 
ergetic fluid originating near the bottom grid and being 
turbulently advected into the observation volume. Simi- 
larly, fluid carried down from the more quiescent region 
above the detection volume has low energy and a smaller 
structure function. 

Figure |9|b) shows the compensated Eulerian struc- 
ture functions, similar to figure |7|^b), but reveals a novel 
insight. Stepping through the vertical velocity bins is 
equivalent to stepping through the energy cascade. Fluid 
coming directly upward from the bottom grid (symbol +) 
carries energy that was recently injected into the large 
scales. As a result, the compensated structure function 
for upward moving fluid is biased towards the large scales. 
Fluid that has downward vertical velocity (symbol A) 
comes from the center region far away from the grid. It 
has had more time to mature, and in this process the en- 
ergy is transported to smaller length scales. Conditional 
structure functions appear to be an effective tool to eval- 
uate whether or not a turbulent flow is fully developed 
and has established a stable cascade. 



4- Third order Eulerian structure functions conditioned on 
the large scale velocity: center region 

Figure [T0| shows the third order structure function in- 
dividually compensated and conditioned on Eu^ in the 
center of the tank. Convergence of third order statis- 
tics was more difficult, so elimination of the two extreme 
conditions was required. The third order structure func- 
tion proves to be similar to the second order in separa- 
tion, symmetry, and collapse to a single functional form. 
The energy dissipation rates found for the three con- 
ditions are: = 25.2cm^/s^, = 21.7cm^/s'^,eA = 
28.7cmVs3. 
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FIG. 9: (Color online) Second order velocity structure func- 
tion conditioned on particle pair vertical velocity (Z direction) 
in the region near the bottom grid. The condition with the 
largest downward velocity has been eliminated due to lack of 
statistical convergence. Symbols re prese nt the following par- 
ticle pair vertical velocities T,Uz/\/ (u'j): + — 3.8 to 2.3, * = 
2.3 to 0.75, o = 0.75 to -0.75, A = -0.75 to -2.3, a) Uncom- 
pensated structure functions b) Individually compensated by 
the energy dissipation rate for each conditional data set. 



F. A Powerful Method for Plotting Conditional 
Structure Functions 

1. Eulerian structure functions conditioned on the large 
scale velocity: center region 

An alternative, and in many ways a more powerful, 
method of visualizing the same data is presented in 
Fig. [TT] Here we show the second order Eulerian struc- 
ture function conditioned on the vertical component of 
the large scale velocity (the same data as Fig. |7|. How- 
ever, the scaled vertical pair velocity is plotted on the 
horizontal axis with conditioned structure functions on 
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FIG. 10: (Color online) Third order velocity structure lunc- 
tion plots conditioned on particle pair vertical velocity and 
individually compensated for each conditional data set. Data 
is taken in the center region of the tank, and the extreme 
vertical velocity plots have been eliminated due to lack of 
statistical convergence . Sym bols represent the following ver- 
tical velocities T^u^/ {ul): * = 2.5 to 0.84, o = 0.84 to to 
-0.84, A = -0.84 to -2.5. 



FIG. 11: (Color Online) Eulerian second order conditional 
structure function versus large scale velocity. Data taken in 
the center region. Each curve represent the following separa- 
tion distances rj-q: + = to 40, * = 40 to 70, o = 70 to 110, 
A = 110 to 140, □ = 300 to 370, x = 370 to 440. 



Eulerian structure functions conditioned on the large 
scale velocity: higher Reynolds number 



the vertical axis. When the structure functions are scaled 
by their value at Sw^ = 0, we find very good collapse of 
the data. The fact that these curves for different r/rj col- 
lapse so well is a striking demonstration that the large 
scales affect all length scales in the same way. The fact 
that the conditional structure functions vary by a factor 
of 2.5 demonstrates the strong dependence on the large 
scales. Note that for Gaussian random fields, the plot in 
Fig. 11 would be flat, and a nearly flat result is observed 
in DNS and grid turbulenc^^. 



Figure 12 shows the effect of increasing the Reynolds 
number. This data is at the center of the tank with the 
grids oscillating at 5 Hz which increases Rx to 380. The 
collapse of the structure function remains. The curvature 
in this flgure is not significantly different from the lower 
Reynolds number data in Fig . 1 1 1 1 indicating that if there 
is a Reynolds number dependence it is weak. 

Figure [l3] shows a comparison of our data with data 
taken in the atmospheric boundary layei^^^ with R\ > 
10*. Atmospheric boundary layer turbulence shows a 
similar collapse of conditional structure functions at all 
length scales. The curvature is also similar in both data 
sets, indicating that the dependence on the large scales 
is similar even at these very large Reynolds numbers. 



In Fig. [TT] it may be expected that the structure func- 
tion at the largest length scales ( x ) are a function of the 
large scale velocity. We see here that the dependence is 
a steep parabola. What is now more clear with this plot- 
ting method is the extent to which all the smaller scales 
are also affected by the large scale velocity; in fact, all 
length scales collapse nearly perfectly onto one parabola. 
The large scale velocity affects all length scales in nearly 
the exact same way, all the way down to the dissipative 
range. 



3. Lagrangian structure functions conditioned on the large 
scale velocity: center region 

Figure [14] shows the Lagrangian structure functions 
plotted versus the large scale velocity, comparable to the 
Eulerian data shown in Fig. [TT] The parabolic shape 
remains, but the curvature is greater for all Lagrangian 
time scales than it is in the Eulerian data. All time scales 
are affected by the large scale velocity. To determine the 
effect of measurement volume bias, we have done this 
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FIG. 12: (Color Online) Eulerian second order conditional 
structure function versus large scale velocity. Data taken in 
the center region at higher grid frequency, 5Hz, resulting in 
higher Taylor Reynolds number 380. Symbols represent the 
following separation distances r/rj: + = to 50, * = 50 to 
100, o = 100 to 150, A = 150 to 200, □ = 310 to 420, x = 
420 to 520. 



FIG. 13: (Color Online) Eulerian second order conditional 
structure function versus large scale velocity. The thin plots 
are from atmospheric boundary layer data'^r/?^: * ~ 100, 
A ~ 400, □ ~ 1000, X ~ 1250. The thick Une is from fig.[TT| 
which has been overlaid for comparison, r/rj: o = 70 to 110. 



analysis for artificially restricted measurement volumes. 
By decreasing the volume by a factor of 2, we observe 
the large r curves to shift by approximately the devia- 
tions between the curves. We conclude that the bias does 
not have a significant effect on the conditional depen- 
dence shown in Fig. [14] for the time differences presented 
(r < 10r„). 



4- Eulenan structure functions conditioned on the large 
scale velocity: near grid region 



Figure [15] shows conditioned Eulerian structure func- 
tion similar to figure [TT] but measured in the inhomoge- 
neous region near the grid (NG) . The structure functions 
here are strikingly different than in the center. The min- 
imum is shifted by more than one standard deviation to 
the left. The inhomogeneity breaks the up-down sym- 
metry so that fluid coming directly up from the bottom 
grid is markedly different then fluid coming down from 
the more quiescent region above (analogous to the * and 
A separation in Fig.]9|. It follows that fluid with an up- 
ward velocity has higher energy than fluid with the same 
velocity magnitude in the downward direction. The at- 
mospheric boundary layer data in Fig. [13] also shows this 
effect with a minimum at Eu^/-^ (w^) = —0.5, presum- 
ably as a result of weaker inhomogeneity. Also notable 



is that the collapse of plots for various r values is not 
as complete as in the central region. This is consistent 
with Fig.[9|b) which shows that the conditional structure 
functions have somewhat different r dependence. 



5. Lagrangian structure functions conditioned on the large 
scale velocity: near grid region 



Figure 16 shows a Lagrangian structure function taken 
in the near grid region, similar to the Eulerian data in 
Fig. [15] The minimum is shifted to the left here also as 
a result of the inhomogeneity in this region of the flow. 
The conditional dependence on the large scale velocity is 
again somewhat larger than in the Eulerian case, and the 
collapse at different time scales is not as complete. 



6. Third order Eulerian structure functions conditioned on 
the large scale velocity: center region 

The third order Eulerian velocity structure function 
plotted versus the large scale vertical velocity is shown in 
Fig. [TT] using data from the center of the tank. Statistical 
convergence is weaker than the second order which limits 
the large scale velocity range available for analysis. The 
collapse seems similar to the second order case shown 



in Fig. 11 although the measurement uncertainties are 
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FIG. 14: (Color Online) Lagrangian second order conditional 
structure function versus large scale vertical velocity. Data 
taken in the center region. Symbols represent the following 
r/r^: + = 0.42 , * = 1.3, o = 3.5, A = 10. 



larger here. The curvature seems to be slightly larger for 
the third order than for the second order case. 
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FIG. 15: (Color Online) Eulerian second order conditional 
structure function versus large scale velocity. Data taken in 
the near grid region of the tank. The structure function is 
heavily influenced by the bottom grid which has skewed the 
symmetry of the plot minima in the negative direction. Sym- 
bols represent the following non-dimensional separation dis- 
tances r/ri: 4- = to 50, * = 50 to 110, o = 110 to 160, A = 
270 to 320, □ = 330 to 450, x = 450 to 560. 



7. Second order Eulerian structure functions conditioned 
on the velocity magnitude: center region 

The second order Eulerian structure function plotted 
versus the magnitude of the pair velocity is shown in 
Fig. [18] using data from the center observation volume. 
The magnitude of the pair velocity is also a useful indi- 
cator of large scale activity. It has no preferred direction, 
and it is a more direct indicator of the instantaneous local 
energy. A similar dependence remains as in Fig. |11[ the 
collapse seems similar, and the curvature is significantly 
larger. 

G. Discussion 

We have provided a comprehensive set of measure- 
ments that shows signatures of the current state of the 
large scales on the inertial range and small scales in tur- 
bulence. Here we wish to discuss factors that might be 
responsible for the dependence of structure functions on 
the instantaneous large scale velocity. First we will ad- 
dress a possible concern that the conditional dependence 
may be a kinematic correlation. Then we will discuss pos- 
sible properties of the large scales that could be impor- 
tant including Reynolds number, anisotropy, mean shear. 



inhomogeneity, and large scale intermittency. We do not 
claim that this list is exhaustive, but it seems that iden- 
tifying the major contributors will be valuable. 

A reasonable suspicion might be that the observed de- 
pendence is a kinematic correlation, meaning that par- 
ticle pairs with large velocity may also have a large ve- 
locity difference simply because the same measurements 
are used in both cases. HosokawaPS identified that Kol- 
mogorov's 4/5ths law requires that velocity sums and 
differences be correlated so that 

«A«_) = I (3) 

where u_ is half the longitudinal velocity difference and 
is half the sum. (For comparison, we have used 
Aur — 2u_ and below Ewy — u+.) KhomyanskjSS! et 
al provide an experimental confirmation of this and in 
a more recent papei^^ present a list of kinematic rela- 
tions. However, several lines of evidence indicate that 
kinematic correlation does not account for the majority 
of the dependence we observe. 

First, note that two independent random samples with 
identical Gaussian distributions have a difference that is 
uncorrelated with the sum, so that the conditional de- 
pendence seen in Fig [11] would be flat. This remains true 
for velocity differences and sums from Gaussian random 
fields. Both of these results can be obtained by consid- 
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FIG. 16: (Color Online) Lagrangian second order conditional 
structure function versus large scale vertical velocity. Data 
taken in the near grid region of the tank. Symbols represent 
the following r/r,,: + = 0.94, * = 2.8, o = 8.0. 



ering the joint pdf of the two samples, and then rotating 
45 degrees to the coordinate system of sums and differ- 
ences. Because the samples are interchangeable, the sum 
and difference axes have to be principle axes of the joint 
gaussian pdf, and the conditional variance of the differ- 
ence is independent of the sum. 

Of course, turbulent velocities are not joint Gaussian. 
However, from the kinematic relations in the literature 
we have not been able to derive predictions for the con- 
ditional structure functions that we consider or for the 
correlation ((Su)^(Au)^) that would capture the main 
conditional dependence we see. 

To make an experimental estimate of the effect of 
kinematic correlations, we conditioned the velocity dif- 
ferences on several other quantities. For each particle 
pair, we calculated the longitudinal and transverse com- 
ponents of the average velocity of the particles, denoted 
Sm|| and Su_l respectively. We then conditioned the 
longitudinal structure functions on the longitudinal and 
transverse pair velocities instead of conditioning on the 
pair vertical velocity. The idea here is that while condi- 
tioning the longitudinal structure functions on the longi- 
tudinal component could have a kinematic corre- 
lation, conditioning on the transverse component 
should have no kinematic correlation. We found that 
conditioning on had a roughly 30% larger effect 
on the structure functions than conditioning on (YiU±). 
Conditioning on Sm^ should have less kinematic correla- 
tion than conditioning on Suy . So more than 70% of the 
effect remains unexplained by kinematic correlation. We 
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FIG. 17: (Color Online) Eulerian third order conditional 
structure function versus large scale vertical velocity in 
the center region. Symbols represent the following non- 
dimensional separation distances r/rj: -|- = to 40, * = 40 to 
70, o = 70 to 110, A = 110 to 140, □ = 220 to 300, x = 300 
to 370. 



conclude that while kinematic correlation may possibly 
make a significant contribution to the conditional depen- 
dence, the majority of the effect comes from the large 
scales. 

An immediate concern when discussing large scale ef- 
fects is if the oscillating grid flow has a Reynolds num- 
ber insufficient for adequate scale separation, and it is 
this which leads to contamination of the small scale stat- 
ics by the large scales. Evidence points to large scale 
dependence not being caused by limited Reynolds num- 
ber. The comparison in Fig. [13] shows that atmospheric 
boundary layer dataP^ with very large Reynolds number 
(i?A > 10^) has nearly the same dependence on the large 
scales as our flow. Increasing the Reynolds number in our 
flow makes very little difference. Additionally, all length 
scales collapse to nearly the same functional form indi- 
cating that limited separation of scales is not the primary 
factor. Taken together, these lead us to the conclusion 
that merely high Reynolds number alone is not enough 
to create small scales that are statistically independent 
of the large scales. 

Our flow is somewhat anisotropic. The ratio of verti- 
cal to horizontal velocity standard deviations is 1.5:1 in 
the center. The effects of large scale anisotropy on the 
small scales has been studied extensivelj^, but it is not a 
major factor in conditional dependence studied here. We 
have analyzed our data by averaging over particle pairs 
with all orientations, so when the structure functions are 
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FIG. 18: (Color Online) Eulerian second order conditional 
structure function versus magnitude of the velocity pair in 
the center region. Symbols represent the following non- 
dimensional separation distances r lr\: + = to 40, * = 40 to 
70, o = 70 to 110, A = 110 to 140, □ = 300 to 370, x = 370 
to 440. 



conditioned on a quantity with no preferred direction like 



the velocity magnitude (Figure 18 1 there should be very 



little contribution from anisotropy. In fact, we find that 
the conditional dependence on velocity magnitude is even 
stronger than the dependence on the vertical velocity 
component. We also observe the conditional dependence 
remains when conditioned on other quantities without 
preferred directions like Suy, and Sit_L. We conclude 
that anisotropy of the large scales is not a significant 
cause of the conditional dependence we observe. 

Sreenivasan and DhruvefSl attribute the strong con- 
ditional dependence of the Eulerian structure functions 
on the large scale velocity to shear in the atmospheric 
boundary layer. In making this argument, they show an 
important piece of information in their figure 6 which 
shows conditioned structure functions in homogeneous 
turbulence from both DNS and wind tunnel grid turbu- 
lence. The conditional statistics in these homogeneous 
and isotropic flows show no apparent dependence on the 
large scale velocity. However, we conclude that shear 
is not the fundamental property responsible in our flow 
since the oscillating grid flow has a much lower shear but 
produces much the same dependence on the large scale 
velocity. The mean velocity gradient normalized with the 
eddy turnover time is 1.2 in the center of our flow and 
we estimate it is in the range of 5 or greater for their 
atmospheric boundary layer data. There must be some 
other properties that exist in the shear flow, but also are 



important in our flow with small shear. 

Our data clearly shows the role that inhomogeneity 
plays in the observed large scale dependence. Our Eu- 
lerian and Lagrangian data near the grid in Figs. |15| 
and [16] show that the structure functions depend greatly 
on the origin of the fluid being swept into the obser- 
vation volume. Fluid coming from energetic regions of 
the tank have larger structure functions than fluid com- 
ing from more quiescent regions. Inhomogeneity is di- 
rectly responsible for the shift of the minimum in Fig. 15 
away from zero vertical velocity. In the center of the 
tank (Fig. Ill, the inhomogeneity is much smaller, but it 
could be responsible for part of the curvature since both 
fluid coming downward and fluid coming upward would 
be coming from more energetic regions symmetrically. 

However, inhomogeneity alone does not account for all 
of the large scale dependence observed. There is also 
a signiflcant contribution from large scale intermittency, 
and it is possible that this is the dominant contribution in 
the center of the tank. Large scale intermittency has been 
difficult to quantify. It can be deflned as any temporal 
fluctuations in the large scales that occur on timescales 
longer than the eddy turnover time, Lju. 

Fernando and DeSilvEp^show large scale intermittency 
can exist in an oscillating grid flow depending on bound- 
ary conditions. We have observed clear signatures of 
large scale intermittency in our flow. Although we use 
their recommended boundary conditions, the velocity 
distribution in the center of the flow is slightly bimodal 
indicative of switching between two flow states. This ef- 
fect is more prominent in preliminary data we took for 
grid spacings of 66 cm and 100 cm than it is in the data 
for 56.2cm presented in this paper. 

Our measurements show a dependence of the condi- 
tional structure functions on the large scale velocity that 
can not be fully attributed to inhomogeneity, and large 
scale intermittency appears to be the most likely cause. 
The clearest evidence for this comes from conditioning 
the structure functions on the horizontal components of 
the large scale velocity, Eu^; and Euj, instead of on the 
vertical component, Em^. The horizontal midplane {x 
and y directions) is much more homogeneous than the 
vertical axis (z direction). Yet, the structure functions 
conditioned on Suj, or Sita, show a large scale depen- 
dence that is only moderately smaller than for the Sw^ 
condition (85% and 72% of the dependence seen in Su^). 
If the inhomogeneous direction shows similar conditional 
dependence on the large scales as the homogeneous di- 
rections show, then it seems that a large part of the 
conditional dependence must come from fluctuations in 
the large scales, and not directly from inhomogeneity. 
Praskovsky et aZ.E^ attribute large scale intermittency as 
a crucial component of the large scale dependence they 
observe. More work is needed to isolate the effects of 
large scale intermittency on small scale statistics in tur- 
bulence. 

We have largely ignored considerations of power law 
scaling which has been a focus of much of the previ- 
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ous work on this subject. Because of the relatively low 
Reynolds number of our experiment, we can not make 
sensitive tests of scaling. However, our data provides a 
plausible picture about how the large scales should af- 
fect power law scaling. If the data in Fig. 11 collapses 
to a single curve, then the dependence of the conditional 
structure functions on r and Uz are separable and the 
large scale dependence will have no effect on the scaling 
exponents of unconditional structure functions. When 
this type of plot does not collapse as in Figs. [15] and [16] 
then the power law scaling will be affected by the large 
scales. 



IV. CONCLUSIONS 

We study a flow between oscillating grids with 3D par- 
ticle tracking and a novel real-time image compression 
system in order to quantify the effects of various proper- 
ties of the non- universal large scales on the incrtial range 
and small scales. 

We measured the mean and variance of the velocity 
as a function of distance from the grids. The oscillating 
grid motion has produced a weak mean flow as well as 
a region near the grid with high velocity variance that 
falls off quickly to a very homogeneous, lower velocity 
variance, region in the center. This profile has been key 
in the determination of the role of inhomogeneity. 

Conditional statistics were employed in order to mea- 
sure the large scale effects. Second order Eulerian veloc- 
ity structure functions were conditioned on the phase of 
the grid, an obvious source for periodic large scale energy 
input. Results show little dependence of the structure 
functions in the center region and surprisingly little even 
near the grid. 

Eulerian and Lagrangian structure functions were also 
conditioned on the instantaneous large scale velocity. A 
large dependence was found in the center, with the Eu- 
lerian structure function increasing by a factor of 2 or 
more when the large scale velocity is large. The depen- 
dence of the Lagrangian structure functions is somewhat 
larger. Conditioned structure functions show that in the 



center of the tank, all length scales are affected in ap- 
proximately the same way. The region near the grid was 
also analyzed and compared with the region in the cen- 
ter. Near the grid, we found a much stronger dependence 
on the instantaneous large scale velocity for both the Eu- 
lerian and Lagrangian structure functions than we found 
in the center. Near the grid, there are clear signatures 
of the effects of large scale inhomogeneity on the small 
scales. Fluid coming up from the energetic region nearer 
the grid has large structure functions, while fluid coming 
down from the quiescent region in the center has much 
smaller structure functions. The functional form of the 
conditional structure functions are also different indicat- 
ing the different histories of the different fluid. These 
measurements provide a clear picture of the way inho- 
mogeneity affects the small scales of turbulence. 

Plotting the conditional structure functions versus the 
large scale velocity provides a powerful method for visu- 
alizing the effects of the large scales on all length scales in 
turbulent flows. We recommend these plots as an effec- 
tive way to compare the effects of the large scales in differ- 
ent experiments. This has been done for grid turbulence 
and homogeneous, isotropic DN^lSl-^^hich show almost no 
dependence of the structure functions on large scale ve- 
locity. Our oscillating grid flow and high Reynolds num- 
ber atmospheric boundary layer turbulence^ show very 
similar dependence. Comparison of conditional structure 
functions in other flows has the potential to clarify the 
effects of the large scales on small scale turbulence and 
to guide the search for universal properties of turbulent 
flows. 
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